Out-of-equilibrium states as statistical equilibria of an effective dynamics 
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We study the formation of coherent structures in a system with long-range interactions where 
particles moving on a circle interact through a repulsive cosine potential. Non equilibrium struc- 
tures are shown to correspond to statistical equilibria of an effective dynamics, which is derived 
using averaging techniques. This simple behavior might be a prototype of others observed in more 
complicated systems with long-range interactions, like two-dimensional incompressible fluids and 
wave-particle interaction in a plasma. 

PACS numbers: 05.20.-y Classical statistical mechanics 05. 45. -a Nonlinear dynamics and nonlinear dynamical 
systems 



The behavior of many complex nonlinear dynami- 
cal systems results in the formation of spatially or- 
dered structures. Examples include two-dimensional 
incompressible weakly dissipative fluids EJ, two and 
three-dimensional magnetohydrodynamics , planet at- 
mospheres ||, electrostatic potential in inhomogeneous 
magnetized plasmas, galaxy cluster formation in self- 
gravitating astrophysical systems Q, vortices in rotat- 
ing Bose-Einstein condensates ||. Similar effects were 
also numerically revealed Q in discrete lattices, in which 
the modulational instability of a linear wave was shown 
to be the first step towards energy localization, followed 
by the nonlinear interaction among breather-like exci- 
tations. Related phenomena have been experimentally 
reported, such as the appearance of oscillons in vibrated 
granular media |Q . 

Besides all analogies, the reasons for this self- 
organization can be indeed quite different, fn systems 
for which a structure originate from any randomly cho- 
sen initial condition, a natural explanation is of a statis- 
tical nature: see for instance Ref. J|] for self-gravitating 
systems or Ref. 0] for two dimensional conservative flu- 
ids, fn other systems, the structures arise and persist 
as a direct consequence of the nonlinear dynamics as for 
instance in Ref. ||. We present in this letter a simple 
and fully tractable model for which both reasons for self- 
organization, statistical mechanics and nonlinear dynam- 
ics, have to be invoked, fndeed, we will show that Gibb- 
sian statistical mechanics can be safely applied only when 
the correct dynamical variables are singled out . 

The model, with long range forces, exhibits long lived 
out of equilibrium structures and related unusual re- 
laxation phenomena jy], like those encountered in self- 
gravitating systems, plasmas or two dimensional fluids. 
Despite the importance of such systems, a clear under- 
standing of their peculiarities is still an open problem. 

We study the repulsive Hamiltonian Mean Field 
(HMF) model |12], that describes the motion of TV par- 



ticles on a circle. Its Lagrangian is 
N a 2 1 N 



This model is an archetype of mean field models for which 
interactions are infinite-range with size dependent cou- 
pling For a special but wide class of initial condi- 
tions, this model jl4|, |l5) has a very interesting dynamics 
since, contrary to statistical mechanics expectations, a 
localized structure (bicluster) appears (see Fig. Q). Its 
initial formation results from a fast oscillation of the 
medium, that nonlinearly drives, on a longer time scale, 
an average motion of the particles in an effective double- 
well potential. The mechanism is qualitatively reminis- 
cent of the parametrically forced pendulum, first ana- 
lyzed by Kapitza the role of the harmonic forcing 
of the pivot is here played, by the self induced fast os- 
cillation of the original dynamics. As for the Kapitza 
pendulum, we average over the fast oscillation and ob- 
tain an effective Hamiltonian for the slow motion. We 
will show that the final bicluster structure is a statistical 
equilibrium of such an effective Hamiltonian. 

Let's consider the evolution of particles for small initial 
energy and velocity dispersion Jf4|. The analysis relies 
on the existence of two time-scales in the system: the 
first one, defined by the linearized dynamics, is intrin- 
sic and corresponds to the inverse of the "plasma fre- 
quency" LU = y/2/2. The second one depends on the en- 
ergy per particle e = H/N via the relation r = et, with 
e = y2e. As a consequence, we use the following ansatz 
di = 6®(t) + efi(t, t). A multi scale analysis shows that 
the correct choice to describe the small and rapid oscilla- 
tions is fi = [V2A + {t) sin(0 + (t)) + o+(r)] cos(0?+^) + 
[V2A-(T)sm(<j)^{t)) + o_(r)] sin((9° + ip) where <j>±(t) 
are fast variables, and d<p±/dt, A±(t), and a±(r) are slow 
variables. The phase tp is defined via the complex num- 
ber Tkf 2 = TV -1 J2f=i e2l6,J ° = |Tkf 2 |e- 2 ^. We introduce 
this expression in Lagrangian and keep terms up to 
order e 2 , averaging over the fast time t, as described e.g. 
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FIG. 1: Bicluster formation : short-time evolution of the par- 
ticle density in grey scale: the darker the grey, the higher the 
density. Starting from an initial condition with all the parti- 
cles evenly distributed on the circle, one observes a very rapid 
concentration of particles, followed by the quasi periodic ap- 
pearence of "chevrons", that shrink as time increases. The 
structure later stabilizes and form two well-defined clusters. 
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FIG. 2: Spatio-temporal evolution of 50 particles with initial 
energy per particle e = 2.5 10~ . The dotted lines present 
the results given by the effective Hamiltonian (^|) whereas the 
solid lines show the evolution of particles initially between 
7r and 2-7T according to the original Lagrangian (fil): except 
for the fast small scales oscillations, they are almost indistin- 
guishable. On the left (between and 7r), the caustics (thick 
lines) given by formula ([?]) reproduce the "chevrons" . 



in Ref. Dropping the e 2 overall factor, the averaged 
Lagrangian reads 



L 



where 
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del 

dr 



(2) 



L+ = — 



A\$ ± u\ - (M° ± ) 2 T 2M° ± a ± uj 2 ± 



M° = eN- 1 £ v =1 e ie i = M?_ 



(3) 



iMi + is the first mo- 



ment associated to the angles and uj± = \M2\)/2 
will turn out to be the fast frequencies (for details see 
Ref. Jl^] ) . (j)± are now cyclic variables and correspond to 
two quasi-conserved quantities, P± — N A\u) 2 _4>±, which 
are adiabatic invariants of the dynamics. Going from 
the Lagrangian to the Hamiltonian formalism, one must 
first remark that, having eliminated the dependence on 
the time derivatives of A±, a±, there is no Legendre 
transform over these variables. In the absence of con- 
jugate variables to the amplitudes A±, the correspond- 
ing Hamilton's equation are simply given, from the least 
action principle, by 8a ± H s i ow = 0. This leads to the 
expression for the amplitudes NA\ = P±to± 3 . Finally, 
from the equations d a± H s i ow = and 4>± = dp ± H s i ow , 
we determine the shifts a± and the phases <p± . 



Rewriting the action, taking into account the relations 
for A± and a±, we end up with a very simple effective 
Hamiltonian that retains only the slow motion 
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1 - |M a 



(4) 



and which describes the evolution of the full system. The 
constant P + and P_ are determined by the initial condi- 
tion. 

For the sake of simplicity, we restrict here to the case 
in which only one mode, P_, is initially excited (P+ = 
0). Dropping the subscript for P, we then consider the 
Hamiltonian 



H eff = H ~y 



P 



1 - | Af a 



(5) 



The particle dynamics determined by this Hamiltonian 
perfectly compares with the one given by Lagrangian (jl]), 
as shown in Fig. 

In order to explain the formation of the bicluster and 
the chevrons observed in Fig. |[ let us first remark that 
the equations of motion of the effective Hamiltonian H e f t 



AV 2 (1-|M 2 |) 



sin 2(0? +V) 



(6) 



are pendulum-like. The evolution of the system thus con- 
sists of N pendulum trajectories in a double-well poten- 
tial, which is self -consistently modified since M2 depends 
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on the 6®. If the initial velocity dispersion is small, sin- 
gularities appear in the density profile which becomes 
infinite along the envelope (or caustics) of the trajecto- 
ries; the chevrons of Fig. [I] are a manifestation of these 
singularities [[TJ)]. This picture allows us also to make 
quantitative predictions. To simplify the calculations, we 
assume that position and depth of the double well poten- 
tial are fixed; this amounts to take ip = const in Eq. (j^) 
(as suggested by Fig. [l]) and \M 2 \ = const which is the 
simplest hypothesis we can make in order to character- 
ize the chevrons. We can compute the time t s of the 
first divergence of the density: it is the time when the 
particles in the bottom of the two wells of the potential, 
that perform a quasi-harmonic motion, cross. We get 
t s = ■Ky/ Nuj-/{%Pe), where P and e depend on the ini- 
tial condition, and lo- depends on \M 2 \. Approximating 
the trajectories by straight lines near the bottom of the 
wells, we obtain the shape of the "chevrons", resorting to 
standard methods |2(J of curve envelope calculation (see 
details in p^| ). The lowest order approximation of the 
n— th shock is 



(t - (2n - l)t s ) 3/2 
V2n - 1 



(7) 



where the \j\pn factor explains the shrinking of the 
chevrons, whereas the 3/2 scaling law is generic accord- 
ing to catastrophe's theory fl9| . Taking for \M 2 \ its mean 
value (0.51 for the present initial condition), Fig. || em- 
phasizes that the agreement with the numerics is quan- 
titatively excellent. 

Since the short time dynamics and the "chevrons" for- 
mation is now clarified, let us next consider the stabi- 
lization of this out of equilibrium state: this is where 
statistical mechanics comes into play. 

Due to the special form of the effective potential, which 
depends only on the global variable \M 2 \, the statisti- 
cal mechanics of the effective Hamiltonian (||) is exactly 
tractable. Indeed, the density of states Sl(E) can be ex- 
pressed in terms of the density of states corresponding 
to the kinetic part of the Hamiltonian Qkin and of the 
angular configurational part Q C onf as follows 

Sl(E) oc [ d\M 2 \n km (E-V(\M 2 \)) Q conf (\M 2 \) (8) 



where V(|M 2 |) = Py/{1 ~ \M 2 \)/2. Using an inverse 
Laplace transform and the Hubbard-Stratonovitch trick, 
one also obtains Q C onf after some calculations. Q(E) in 
Eq. (SJ) can thus be evaluated Q by the saddle point 
method, which leads to the equilibrium value \M 2 \ as 
a function of E/P. This explains the numerical result 
\M 2 \ ~ 0.5 found by the authors of Rcf. Q (their initial 
conditions correspond to E/P = y/2/2 and \M 2 \ = 0.51). 
In addition, this shows that other initial conditions lead 
to other values of \M 2 \, which are in excellent agreement 
with numerical simulations (see Fig. ||). The long life- 
time of these out-of-equilibrium states is therefore fully 
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FIG. 3: The solid curve shows the theoretical prediction of 
\M% | as a function of the ratio between the energy E and the 
adiabatic invariant P, whereas the filled circles correspond to 
numerical results for the full Lagrangian (hi). 



understood, since they appear as equilibrium states of an 
effective Hamiltonian that well represents the long-time 
motion. 

Once |MJ| is known, it is easy to show that the equi- 
librium velocity distribution is Maxwellian with temper- 
ature T = 2(E C )/N, and that the distribution of angles 
has a Gibbsian shape p{6) oc e~ v ^l T , with the potential 
inferred from the equations of motion. However, whereas 
\M 2 \ quickly reaches its equilibrium value, the distribu- 
tion p(9) relaxes very slowly, in the numerical experi- 
ments. Moreover, the relaxation time increases with N . 
Nevertheless, the density profile p{9) obtained in long- 
time simulations with the full system ([!]) fully agrees 
with the one obtained with the effective dynamics (JEI), 
as shown in Fig. |I| This makes this model a good candi- 
date to study the unusual relaxation properties and the 
non equilibrium states observed in other systems with 
long range interactions [ pT| . 

We thus conclude that the dynamics of the effective 
Hamiltonian <S) parameterizes very well the one of the 
full Lagrangian (|]), for short as well as for long time. 
The variational multi-scale analysis we used allows us to 
exhibit naturally the adiabatic invariants and to preserve 
the Hamiltonian structure of the problem (to leading or- 
der), making it well suited for a statistical treatment. 
This leads us to predict statistical properties of the full 
system, as the asymptotic value of \M 2 \. Moreover, the 
effective Hamiltonian gives the opportunity to study nu- 
merically the relaxation towards equilibrium of the bi- 
cluster, whereas it was not possible in the real dynamics, 
because of computational limitations. Indeed, let us ob- 
serve that the ratio of the typical time scale of the two 
dynamics is of order 100 or even larger at smaller ener- 
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localization in classical and quantum mechanics. 
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FIG. 4: Comparison of the non Gibbsian density distributions 
p(0) of particles obtained with the original Lagrangian ([j]) 
(circles) and with the effective one (pi) (solid line). Both re- 
sults have been obtained for N = 10 particles and are aver- 



aged over times corresponding to r = 10 3 
per particle is e = 2.5 10~ 5 . 
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gies. 

We have thus described a simple mechanism to explain 
the existence of a stable out-of-equilibrium structure in 
a Hamiltonian mean-field model. This model deserves 
special attention for different reasons. It is probably 
the simplest iV-particles solvable model in one dimension 
which exhibits such a stabilization effect, corresponding 
to the coupling of very fast oscillations self interacting 
with a slower motion. — This model is moreover a sim- 
ple analogue of other examples of nonlinear interactions 
of rapid oscillations with a slower global motion like the 
piston problem |2]| : averaging technics could be applied 
to the fast motion of gas particles in a piston which itself 
has a slow motion p2fl . Examples can also be found in 
applied physics as for instance wave-particles interaction 
in plasma physics |2j| , or the interaction of fast inertia 
gravity waves with the vortical motion for the rotating 
Shallow Water model 
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